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Abstract 



Long-established results for the low-energy photon-photon scattering, 77 — >■ 77, have recently 
been questioned. We analyze that claim and demonstrate that it is inconsistent with experience. 
We demonstrate that the mistake originates from an erroneous manipulation of divergent integrals 
and discuss the connection with another recent claim about the Higgs decay into two photons. We 
show a simple way of correctly computing the low-energy 77 scattering. 

PACS Numbers: 13.60.Fz, 13.40.-f, 12.15.Lk 

1 Introduction 

After Dirac proposed the theory of negative energy solutions of his equation [IJ, it was realized that 
photons can interact with other photons by polarizing the vacuum. Photon-photon scattering was 
qualitatively considered in this context by Halpern [5], and its cross section, for the case of photon 
energies low compared to the electron mass, was determined by Euler and Kockel in 1935 [3111]. If the 
energy of each of the colliding photons is lo in the frame in which their total momentum vanishes, the 
low-energy differential cross section is 



where a ~ 1/137 is the fine structure constant and m is the electron mass. 

High energy scattering was considered soon afterward [51 IB]- A thorough analysis of the scattering at 
all energies, including partial cross sections for various polarization states, was carried out in [71[S], 
using the then new diagrammatic technique of Feynman. Since then, the photon-photon scattering 
cross section has been confirmed with other methods, and even higher-order QED corrections have 
been computed |9]- Results obtained up to 1971 are reviewed in [TU] and more recent developments 
are summarized in |11| . 

Very recently, the classic result for the low-energy cross section ^ has been questioned |121 [T3] . In 
those papers, the cross section is found to be many orders of magnitude larger, since it is not suppressed 
by powers of (w/m), but is proportional to 1/w^, 



As we will demonstrate in this paper, this claim is incorrect. It has already been pointed out [14! that 
it contradicts existing laboratory bounds on the the photon-photon cross section, obtained by colliding 
laser beams. We show in addition that a cross section increasing with the inverse squared energy 
of the colliding photons limits the mean free path of visible light due to collisions with the cosmic 
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microwave background radiation (CMBR) to less than the distance between Earth and Jupiter. Thus 
the fact that we can sharply see much more distant astronomical objects proves that the low-energy 
photon-photon scattering must be significantly suppressed, as predicted by eq. ([Ij. 

The matrix element for the photon-photon scattering is absent at the tree level since photons are 
neutral. It arises only at the loop level. The sum of all contributing loop diagrams must be finite 
since there is no parameter in the QED Lagrangian whose renormalization could absorb a divergence. 
Refs. [T^ [13] found an incorrect result because of the assumption that if the sum of those diagrams is 
finite, they can be calculated without any regularization. In fact, even though the sum of the diagrams 
is finite, each of them separately is divergent. Calculating the sum is somewhat delicate and is easiest 
done with regularized loop integrals (see, however, an alternative calculation in |15| and another point 
of view on avoiding regularization in |16|). 

Interestingly, a similar error pTl [T5] has recently cast doubt over the rate of the Higgs boson decay 
into two photons. That process, too, is loop induced, and the sum of contributing loops is finite. 
But individual loop integrals are divergent and must be regularized, as has already been thoroughly 
discussed in this context [T9 l [20 l [2T | [22 l [23] . 



2 Mean free path of photons in a microwave background 

The CMBR is a gas of photons with the spectrum of a black body at a temperature of 1//3 = 2.725 
K. Here we want to compute how far a visible- light photon with energy ~ 2.5 eV can travel in 
such a gas before scattering, from the point of view of an observer in whose frame the CMBR is 
isotropic. (We will call it the LAB frame. For the purposes of this discussion an Earth-based observer 
is a good approximation.) Consider one mode of the CMBR radiation, characterized by its energy 
E and inclination angle 9 with respect to the direction from which the visible photon is incident. 
The relative velocity of the two photons (as seen in the LAB frame) is vi + V2 = {1 + cos 9, — sin9), 
Wi + W2I = 2 cos I (we use the units c = h = ks = 1)- In the frame where the total momentum of the 
photons vanishes, each has the energy lo given by 

1 / 9 

uj = -■^2EE^ (1 + cos 9) = cos- y/EE^. (3) 

That energy determines the scattering cross section. Collisions with photons in this particular mode 
will occur at the rate 

dTEe^\vi+V2\<jdp{E) (4) 

where 

E^ di^dcosg 
'^"^'^)=2^ exp(/3g)-l 

is the density of CMBR photons with energy E, and a is the scattering cross section. Integrating over 
the energies and directions of the CMBR photons we find the mean free path. Between collisions, the 
visible-light photon will travel on average the distance 



/I Q j^2^ 
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We now consider the two formulas for the low-energy cross section. If we use the classical result ([T]), 
we find the total cross section 
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and the mean free path 
A = 
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Using m = 0.511 MeV we find A ~ 7 • 10^* meters, a distance that would take hght about 10^^ times 
more time to travel than the age of the Universe. In other words, the CMBR is a rather transparent 
medium at visible frequencies. 

However, if we take instead the cross section suggested in [T^ I13| . we find from eq. ^ 
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which gives a much shorter mean free path, 
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or ApK — 3 • 10^^ meters, equivalent to about 15 light minutes. For comparison, the orbital radius of 
Jupiter is about 8 • 10^^ meters. If the mean free path of the visible light were so much shorter than 
even the radius of Jupiter's orbit, no stars would be visible on the night sky. Clearly, the result eq. ^ 
is at odds with experience. 

The situation with eq. ^ is actually even worse. Since the cross section is not suppressed by the mass 
of the electron, there would be additional positive contributions from other charged fermions that 
would differ only by the coupling constant, and would further decrease the mean free path. This lack 
of suppression by the inverse mass of the loop particle contradicts the Appelquist-Carazzone decoupling 
theorem [24]. 



3 Determination of the photon-photon scattering 

In this section we present a derivation of the photon-photon scattering matrix element in two reg- 
ularization schemes: dimensional and Pauli-Villars. We consider the box diagram shown in Fig. [TJ 
External photons carry momenta /ci, . . . , fc4 which we will consider as incoming, /ci + fc2 + ^3 + ^4 = 0. 

There are six ways in which the four momenta can be arranged around the oriented electron loop. 
However, diagrams that differ only by the direction of the electron line give identical results so it is 
enough to compute three of them, corresponding to three cyclic permutations of /ci , fc2, and k^. (If there 
was an odd number of photons coupling to the electron loop, the diagrams differing by the direction 
of the electron would cancel one other, resulting in a vanishing amplitude. This is the theorem due to 
Furry [25J.) 

At low energies of external photons it is especially easy to compute the diagrams in Fig. [TJ We 
simply Taylor-expand the electron propagators in the external momenta, so that each propagator's 
denominator becomes simply (^ — to) ^ = {i + m) / (q^ — to^) where q is the loop momentum. Such 
expansion does not lead to any spurious divergences, and commutes with the integration over q. We 
now explain how this integration is performed in two regularization schemes. 
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Figure 1: Virtual electron loop inducing the four-photon coupling. 
3.1 Dimensional regularization 

Now that q is present in the denominators only through , also in the numerator we can replace all 
scalar products of q with other vectors by powers of q^ times products not involving g, 



g^i . . . g^^" 
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Here D is the space-time dimension and S ((ji'^i'^^ _ g,At2n^iA'2n ~) jg Qf products of n metric tensors 

g, totally symmetric in all indices /i^; it has (2n — 1)!! terms. Terms odd in q vanish upon integration. 

The powers of q^ resulting from (llip can be canceled against the denominators and the loop integration 
can be completed using 



A^q 



{2T: f (q^ ~rn^ + W (4^) 



D/2 



r(«) 



(12) 



Each of the three diagrams contains terms with the exponent a = 2, leading to a divergence F (2 — -j) ~ 
1/ (_D — 4). The divergences cancel when we add all three contributions. But individual diagrams 
containing singularities 1/(1? — 4) have also Z?-dependent factors, arising from the averaging in eq. 
(jlip . The resulting finite contributions do not cancel among themselves. 



How do these remaining terms depend on ml We remember that they arise from the a = 2 sector, 
therefore they scale like mP (the overall dimension of the 77 — !■ 77 amplitude) . There are other terms 
that scale with this power, arising from convergent integrals like J d^q/ (q^ — m^) . The essential 
point is that the sum of all mP terms, including the remnants of singularities, adds up to zero. The 
total result for the amplitude turns out to be suppressed by four powers of 1/m. 



3.2 Pauli-Villars regularization 



Another way of carrying out this calculation is to stay in four dimensions but add another amplitude, 
with the electron replaced by a very heavy particle of mass M, and with an opposite sign than the 
electron loop. The calculation proceeds very similarly to the case of dimensional regularization, with 



two changes. In averaging over the loop momentum directions (|lip we replace 



r(g+2n-l) 

r(-D) 



by its value 



at D — A, {2n + 2)!/6. The formula for the loop integration ([12]) is also replaced by its _D = 4 value, 
except in the divergent case a = 2. In the dimensional regularization, this divergent integral gives 
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In TO . In the Pauli-Villars approach one finds a convergent combination 
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When the three diagrams are added, this logarithm cancels, but now there are no finite remnants of 
the singularities. Instead, the m° terms from the convergent diagrams are canceled by the terms 
from the Pauli-Villars subtraction. Since they are independent of the electron mass, they are the same 
in the amplitudes with the electron and with the very heavy particle, and cancel in the difference. 

In both regularization schemes, the only remaining result is suppressed by the electron mass. 



3.3 Potential error from neglecting regularization 

We have just seen that the regularization is crucial in computing the photon-photon scattering ampli- 
tude, even though the final result does not contain divergences. We now want to inspect more closely 
the part of the amplitude that does not contain external photon momenta, and thus scales like the 
zeroth power of the electron mass. 
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with 
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gl^upa ^ gf'^qPq" + five other terms. 



where the terms not shown in 5*2 have the other five distributions of indices so that both Si and 5*2 
are totally symmetric in fi, p, a. The second line in (|14[) contains four powers of the loop momentum 
q and thus represents divergent integrals. Without regularization, these divergent integralssimply do 
not have a meaning. If we apply the averaging procedure to these terms, we find 

{q''q''q''qn = ■^^TYf''^'^ ' ^^^^ 

so that if D = 4, the second line of vanishes, as does the term ^ m? in its first line. Thus, if 
the regularization is neglected, one is left with only the first term m^Si which, after the q integration, 
gives a result independent of the electron mass, scaling like wP , 

iMm° = -^a^-Sf ''°'ei^e2^e3pe4cr (16) 

= -^a^ (ei • £2 £3 • £4 + ei • 6362 • £4 + ei • e4e2 • £3) , (17) 

where are the polarization vectors of the four photons. This dependence of the amplitude only on 
the polarization vectors (and not on the photon momenta) means that the induced coupling of the 
photons involves only their vector potential (the induced effective operator is proportional to (A^)^), 
and not its derivatives. It is not possible to construct such a coupling in a gauge invariant way. 

This violation of gauge invariance may also generate photon's mass. For example, if two of the external 
photon lines in Fig. [1] are contracted, the resulting two- loop diagram generates an operator ^ A^, thus 
giving the photon a mass. 
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In order to see how the cross section in eq. ([2]) follows from the amplitude (fT7|) . we define two transverse 
polarization vectors for each photon, e^'^, with e* perpendicular to the scattering plane and lying 
in that plane. We do not include here the longitudinal photon polarizations, present if the photon 
becomes massive, even though they may dominate the cross section; however, our goal here is merely to 
explain how the result ([2]) is related to the gauge- invariance violating amplitude (|17p . For the eq. (flTl) 
to give a non-zero result, each polarization must be represented an even number of times (otherwise 
there will always be a factor in every term). There are eight possible such polarization assignments, 
giving the following values of the three terms in (1171) , 

Mini ^ 1-f 1-1-1, 

7^2222 - 1 + COS^ 6* -f COS^ 0, 

A^1122 - I +0 + 0, 

A^i2i2 ^ -f COS 6* -I- 0, 

A^i22i ^ + 0-FCOS6', 

and the last three amplitudes enter with a weight factor of 2, due to the symmetry 1 -H- 2. The various 
amplitudes differ by the polarization of some photons, so they do not interfere. The sum of their 
squares gives 3^ -|- (l -I- 2 cos'^ 9) +2 + 4: cos^ 9 = 4(3 + 2 cos^ 9 + cos^ 9) , the angular structure of the 
(incorrect) result quoted in ([2]). In fact, a sum over the polarizations of the final state photons and an 
average over the polarizations of the initial state photons, leads to the cross section given in eq. 

What went wrong in the above procedure? The formulas (jlSp cannot be applied to the divergent 
integrals in the second line of (fT4|) in = 4, without regularization. If we stay in D dimensions, the 
terms we found to be zero in the — 4 limit given finite contributions that cancel against the first 
term of the integrand, '-^ m'^Si. In this correct treatment the resulting amplitude is suppressed by four 
powers of 1/m. 

The recent incorrect claim about the decay H — > 77 |17[ 118) originated with a similar, but somewhat 
simpler integral. An example of a contribution to that process is shown in Fig. ([2]). There are only 
three propagators, and the divergent integrals are present in the combination |18] 



V (D) = J d^q 



[q^ — ni^ + iO)^ 



Without dimensional regularization, if we take _D = 4, it seems that this integral vanishes after 
averaging over q with help of (jlip . As we have seen with the example of 77 scattering, and as has 
already been discussed in the literature [2Tl[T9l[T^, such manipulations with unregulated, divergent 
integrals are unjustified. In the case of the Higgs decay, they lead to the incorrect conclusion that 
I^u (D 4) vanishes. In fact, in the limit of a very heavy Higgs boson, the correct finite result of /^i, 
gives the most important contribution. 




Figure 2: An example of aW boson loop mediating the Higgs boson decay into two photons. 
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3.4 Results for polarized photons 



The correct result of the loop integration in the 77 — > 77 amplitude contains scalar products among 
photon momenta ki, in addition to the polarization vectors e^. The effective photon-photon coupling 
induced in this way is described by operators involving the electromagnetic field tensor and is gauge 
invariant. We know calculate the scattering cross sections for various polarization situations. Instead 
of the linear polarizations we have just considered, the scattering amplitudes will be presented in terms 
of circular polarization states. Thus we introduce 

describing right- and left-handed polarization states, respectively. There are four possible initial polar- 
ization states, but it is sufficient to consider just two of them, ++ and H — . We get three independent 

scattering amplitudes, , and We describe kinematics in terms of Mandel- 

stam variables s = (A:i -|- ^2)^ = 4:E'^ and t = {ki + k^)"^ = —2E'^ (1 — cos9) and find 
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(18) 



We note that the amplitude Ai++^ vanishes at the leading order in E/m expansion at which the 

other amplitudes are finite. In order to compute it, we have to evaluate two more terms in the Taylor 
expansion. 

All other amplitudes can be obtained from eq. (|I8p using space and/or time reversal and the crossing 
symmetry. For example, iM+-+- = - ""!f and iM+ — + = -^^'^^fi^ 



3.5 Total cross section 

Once the polarized amplitudes have been evaluated, the unpolarized cross section can be easily found. 

We quote here only the leading low-energy result (thus we neglect A^+4__| and seven amplitudes 

related to it) for the cross section averaged over initial and summed over final polarizations, 

da (77 ^77) ^ 1 \M++++\' + \M++-\^ + \M+-+-\' + \M+—+\' 

dO. 2567r2w2 2 ^ ' 

(3-|-cos^6l) , (20) 



(I807r)^m8 



in agreement with the classic result ([T]) . When integrated over both and (/> from to tt (we integrate 
only over one hemisphere since the two final-state photons are identical), this gives the total photon- 
photon scattering cross section, 

o- (77 ^ 77 = TTTPTT 21) 

in agreement with [53]. Other texts seem to have misprints in these results |27l 128] . 
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4 Conclusions 



We have demonstrated that the recently claimed result for the 77 — > 77 scattering cross-section © 
must be wrong. The photon-photon coupling is induced by virtual loops with charged particles and 
is suppressed at low photon energy by the inverse power of the electron mass. The result ([2]) lacks 
this suppression and yields a very large cross-section, therefore a short mean free path of visible 
photons even in the rare cosmic microwave radiation background. Such short path would obscure all 
astronomical objects as close as Jupiter. 

We have showed that the error resulted from manipulating unregulated divergent integrals. A similar 
error misled the authors of [171 [TS] in the context of the Higgs decay to two photons. Both processes 
arise only at the loop level and their amplitudes must be finite, since there are no parameters in the 
Lagrangian that could absorb a divergence. However, both processes are usually computed from a sum 
of several diagrams, among which some are divergent. For this reason, a regularization of individual 
contributions is necessary. 

Interestingly, the mistakes in these recent studies of 77 — 77 and ii" — >■ 77 led to confusions about var- 
ious types of decoupling. The correct result for the former process does respect Appelquist-Carazzone 
decoupling theorem in the limit of low photon energies or large electron mass, whereas the incorrect 
result of [121 [T^ does not. On the other hand, the correct result for the Higgs decay does not vanish, 
as one could naively expect, in the limit of large Higgs mass (or, equivalently, low W boson mass; 
this type of decoupling affects for example quarks but not the longitudinal W components) , while part 
of the reason why [ITl [18] believed their result was that it did vanish in that limit. 

We have also showed how the 77 — >■ 77 amplitude can be calculated in the low-energy regime, and 
how an expansion in powers of the photon energy to the electron mass ratio can be organized. This 
tutorial illustrates useful techniques of loop calculations: averaging over loop momentum direction, 
loop momentum integration, and various regularizations. We hope it will be helpful for other similar 
loop calculations. 
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